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■ Abstract 

\ We take a closer look at the recently discussed models of hadrons based on holo- 

qq ' graphic ideas and chiral symmetry breaking in five dimensions. We study the two point 

OS ■ correlator in detail and look at the field theoretic properties needed to be satisfied by 

the correlator. It is shown that the spectral density becomes negative, violating the 
basic requirements for the Kallen-Lehmann representation. We briefly discuss the im- 
plications of this violation of the positivity of the spectral density and also discuss 
possible origin of such a violation in these models. Put simply this means that such 
models are not physical descriptions of the hadron spectrum. 
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^ ■ The correct and complete theory describing the confinement of quarks and gluons into 

hadrons, and explaining the hadronic properties from first principles was, and still is, one 
of the major challenges. Though there exist many phenomenological models which seem to 
describe the low energy data reasonably well, a complete theoretical understanding is still 
lacking. The holographic principle or the so called AdS / CFT correspondence pQ has opened 
a new window to look into the problem from a completely new and different perspective. 
Very loosely speaking, the correspondence is a conjecture of the equivalence of the generating 
functional of (large N limit of) certain conformal field theory (CFT) in <i-dimensions and 
the (d + l)-dimensional supergravity (string theory in a particular limit to be more precise, 
though we will use the word supergravity throughout the note) effective action evaluated 
at the boundary. The boundary value of the bulk fields is supposed to play the role of the 
source in the CFT generating functional. If we generically denote the (d + l)-dimensional 
fields as 0, then the correspondence can be schematically summarized as 

ZAdslM = I V(j> exp{-I[(j>)) (1) 
= Z C ft[J = M = (exp{ I d d xO^)) 
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where 0o is the boundary value of the field </> which acts as a source corresponding to the 
operator O in the CFT generating functional. For the (d + l)-dimensional AdS space, the 
line element is 



where x M are the the <i-dimensional coordinates and z parameterizes the extra direction, 
with z = describing the ultra-violet (UV) boundary and L denotes the radius of curvature, 
which we set to unity. 

What it means in practical terms is that a perturbative calculation on the supergravity 
side can be translated into a non-perturbative result in the gauge theory sector and vice- 
versa. If this be so, it simply implies that the seemingly impossible task of getting the 
genuine non-perturbative results and calculating the non-perturbative quantities within a 
gauge theory is no more impossible. Instead, one can approach the problem from the dual side 
- where perturbative calculations, though hard, are possible. This has motivated considerable 
interest in getting the hadron spectrum and properties. Till very recently, most of the 
models constructed and studied can be broadly thought to fall under the top-down category 
of models i.e. the approach was to start from some string (or supergravity) theory and try 
to obtain the low energy description by demanding/imposing certain consistency conditions. 
For some of the earlier works, see [2|. Encouraged by these explorations, non-supersymmetric 
holographic dual models for hadrons have recently been proposed j^j, jlj. Compared to 
previous studies, these models are phenomenological and the approach adopted is bottom-up. 
Guided by the basic ideas of the correspondence principle, one identifies the corresponding 
conserved currents which appear in the gauge theory as a dual description to the fields in the 
supergravity theory. This thus determines the (minimal) field content of the five dimensional 
theory of gravity, as dictated directly by the low energy sector. Thus the name, bottom- 
up. The models are still at the stage of being called toy models and they employ a very 
small sub-set of the possible field content. One can then go ahead and calculate various 
n-point Green functions (or correlators) of the fields in the 5-dimensional gravity theory. 
The correspondence relates such a calculation to the correlators of conserved currents in a 
suitable gauge theory. Using only a minimal sub-set of fields, the authors have shown that 
the models are quite robust and predictive to within 10% accuracy. To capture the essentials 
of chiral symmetry breaking, some specific boundary conditions are imposed on the fields 
and their derivatives on the so called "infra-red (IR) boundary" (zm » zyy). The theory 
is conformal only close to the UV boundary. As one moves away from this UV boundary, 
the theory is no more conformal. In a complete microscopic description, this should be 
incorporated by appropriately modifying the geometry which was AdS to start with. In 
the absence of a complete description, this is done by putting certain artificial and ad-hoc 
boundary conditions and one hopes to capture the essential and broad features. 

In this Letter we take a closer look at these models and explore them in more detail. Given 
the robustness of the predictions within these models, it is tempting to investigate how 
far can they go in describing the experimental observations. The two point functions are 
the simplest and most straightforward to be evaluated within the 5D gravity theory. Their 
importance is not just being the simplest objects calculable, but lies in the fact that they offer 
a dual description to the current- current two point correlator ,(j^j u ), in the gauge theory. 
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This correlator shows up in e + e~ annihilating into hadrons, r hadronic decays, hadronic 
contribution to (g — 2) M and many other places. Moreover, the current-current correlator 
is used to extract the masses and decay constants of the mesons in the large N c limit of 
QCD. Furthermore, it is known that the masses and decay constants obey certain sum rules, 
namely the Weinberg sum rules and generalizations [5J. It is thus important to have a precise 
and accurate theoretical description of the two point current correlation function. 

Before proceeding to explore the two point function in more detail, let us briefly mention the 
field content and other details of the toy model |3] that we will be working with. We follow 
[3] for the notation and general setup of the model. The field content of the 5D theory is 
dictated by the operators in the 4D QCD. In principle, there should be an infinite number 
of fields corresponding to an infinite number of gauge invariant operators in QCD. However, 
for the purpose of chiral symmetry breaking and its essential consequences, it suffices to look 
at a very small number of operators and therefore a small number of fields in the 5D theory. 
As is known, the chiral dynamics is quite effectively described by an SU(Nj) L x SU(Nf) R 
theory. We therefore have the following three operators that are most crucial to this effect 
(corresponding 5D fields are written in front of each of them): 

q L (x)YT a q L (x)^iA^(x,z) 

q R (x)rT a q R (x)^A%(x,z) (3) 
qi(x)q R (x)^(2/z)X«(x,z) 

where T a are the group generators and i, j are the flavour indices. Restricting to two flavours 
implies that i,j = 1,2 and the T a, s are the three Pauli matrices. We are thus looking at a 
chiral SU (2)l x SU {2)r theory. With this minimal field content and ignoring any interactions 
for the time being, the 5D action is 



5(5) = J d^xdz^fg Tr 
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l -(F 2 L + F 2 R ) + \DX\ 2 -m 2 X 2 



(4) 



where F^ R denotes the field strengths for the left and the right gauge fields and Da denotes 
the covariant derivative. 

At the IR boundary, some boundary conditions are imposed on the fields. Moreover, the 
reason behind introducing the IR boundary in the theory is to parameterize the effect of 
chiral symmetry breaking in an effective fashion. Else, one would be forced to examine how 
the geometry of the bulk changes away from the UV boundary and what should be the 
metric far away that will suitably describe chiral symmetry breaking in the 4D gauge theory. 
The introduction of the IR boundary and suitable boundary terms in a phenomenological 
way takes care of this aspect. For the gauge fields, we impose the boundary conditions: 
F£ Z R = and choose to work in the gauge A z = 0. For the field X, the expectation value at 
the UV boundary is the quark mass matrix and the quark condensate will effectively fix the 
other constant (at the IR boundary) in the solution to the equation of motion for field X. 
In principle, the UV boundary corresponds to z = 0. However, in practice, the boundary 
conditions are specified at z = z uv and finally the limit z uv — > is taken. 



3 



Define the vector and axial- vector gauge fields as appropriate linear combinations of Al and 
A R . Let us focus on the vector gauge field Va = a + (Ar) a) and let Va((?, z) denote the 

Fourier transform of the vector field with respect to the 4D coordinates. We have suppressed 
the group index for convenience. In the V z — gauge, the linearized equation of motion for 
the transverse part reads 

d s (U z %(q,z)} + £%(q,z) = U (5) 

The solution is a linear combination of the Bessel functions and can be written as 

%{q,z) = C,{q)qz[bJ 1 {qz)+y 1 (qz)} (6) 

The boundary condition at the IR boundary, namely d z V fl \ z=ZlR = 0, fixes the relative 
constant between the two terms to be 

u _ y^qzin) ^ 



We can now use the correspondence principle to interpret the above solution at z = zjjy as 
the source for a current. For the two point function we get 0, 

VMoM) = U» - ^f) i%) (8) 



Q- 



where 



n(g)| 



z=z vv 




Q 



zuv 



y^qzuv) + bj (qz uv ) 
yi{qzuv) + bJi(qzuv) 



(9) 



This expression can be expanded about zjjv = and we retain the leading terms, which 
leads to 



n(g)| 



z=z uv 




71 

\og(qz uv /2) +7B + -6 



(10) 



The logarithmic term may be identified as the contribution to the current-current correlator 
arising from the lowest order quark bubble diagram. Comparing this expression with the 
one obtained within QCD with iV c colours fixes the 5D gauge coupling, g 5 , in terms of iV c . 

The next task is to estimate the hadron masses and decay constants. This is easily done by 
comparing the two point function obtained above with the corresponding expression that one 
obtains in the large N c chiral theories. In such a theory, the two point function is expressed 
as a sum over narrow resonances. The poles of the two point function yield hadron masses 
and the residues at each pole are the decay constants. The correlator thus has the following 
form (the ie is implicit) 

ntaWE^ (ii) 

n H n 

In particular, the above form and the fact that we have approximated the correlator as a sum 
over narrow resonances imply that the spectral density is a comb of delta functions peaking 
at the hadron masses. The spectral density is nothing but the imaginary part of n(g 2 ). Also, 
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let us recall that the Kallen-Lehmann representation for the two point function in a generic 
field theory implies that the spectral density is a positive quantity [7j. This is an important 
property that it must satisfy and we'll see that this property plays an important role in our 
analysis. 

To estimate the hadron masses, we look for the poles of the two point function. We can 
use either Eq(JHJ) or Eq(|TU|) for this purpose. The poles are given by the zeros of Jo(q z ir)- 
Then using the experimental value for m p as an input fixes Zir. The decay constants are 
the residues at these poles. 

The estimated numbers for the masses and decay constants [3] compare well with the ex- 
perimental values and we get the impression that the model is very robust and quite close 
to reality. However, as we'll see below, this is not the real and complete picture. We take a 
closer look at the two point function. In particular, we are interested in the corresponding 
spectral density, p(q 2 ). Both the vector and axial- vector spectral densities are measured to 
very good accuracy, for example, in the hadronic r decays and are directly related to the 
decay width by optical theorem 0. The data clearly shows the p and a\ resonance peaks and 
supports the theoretical expectation of the spectral densities to be positive functions. We 
use Eq([ll)|l (or equivalently EqtjHJ)) to extract the imaginary part, and therefore the spectral 
density. Using the usual ie prescription, it is easy to obtain the imaginary part which is 
nothing but the following 

yo{qzm) 



p{q) = Im 



Jo((q + ie)z IR ) 



(12) 



We would like to study the above spectral density in more detail and convince ourselves 
that it satisfies all the field theoretic properties, like the positivity condition, before we can 
proceed further and calculate masses and decay constants from the two point correlation 
function. To this end, we consider the following function 

/(*) = y °} X ) = Re[f(x)] + ilm[f{x)\ (13) 
Uoyps i lit) 

Clearly, the imaginary part of the function is nothing but the spectral density itself rewrit- 
ten in terms of a different variable. In Figure 1, we have plotted the imaginary part of the 
function, Im[f{x)\. 1 A quick look at the plot gives the impression that the spectral den- 
sity is indeed a comb of delta functions, as expected and desired. The zeros of the Bessel 
function Jq{x) are the positions of the resonance masses and the residues at these values 
will correspond to the decay constants of the mesons. Let us now take a more closer look 
at the function itself. The function under investigation is a ratio of two Bessel functions. 
Further, the Bessel functions are known to have an oscillatory behaviour, with Jo(x) and 
3V X ) having opposite behaviour with respect to each other. Figure 2 is a plot of the two 
Bessel functions, clearly showing these features. 



1 A quick way to see that the imaginary part will have a profile very similar to a sum of delta functions 

2 4 6 

is to look at the series expansion, Jq(x) — 1 — ^- + |j — 2§o4 + The ie prescription can now be used 

to obtain the imaginary part. For plotting the graph, we have chosen e = 10~ 7 and rescaled the y-axis. It 
should be borne in mind that any other small value of e is equally good and rescaling simply helps in having 
an enlarged picture and does not change the shape and nature of the curve. 
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The relevant quantity (the imaginary part of the function f(x) which is the spectral density) 
is a ratio of these two different Bessel functions and because of the features mentioned 
above, 3V X ) can cross the real axis between two zeros of Jo(x), thereby yielding negative 
values for the imaginary part of the function, and therefore the spectral density. To further 
substantiate our claim, we look at the behaviour of the imaginary part of the function f(x) 
more closely. In Figure 3, we plot Im[f(x)] for various smaller intervals of x and show that 
it indeed acquires negative values. The reason that this feature is not evident in Figure 1 
is due to the very large values that the function acquires close to the resonances. However, 
when we look at the behaviour in regions slightly away from the resonances, the negative 
values show up, as shown in Figure 3. 
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Figure 1: Coarse-grain view of Im[f(x)] as a function of the argument, x 
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Figure 2: The two Bessel functions: J$(x) (left) and y$(x) (right) 

Let us put all the individual pieces of information together to get a final and complete 
picture. We have seen above that a cursory look at the imaginary part gives the impression 
that it is a comb or sum of delta functions - the position of the peaks will give the masses of 
the mesons and the residues at these values will be the respective decay constants. It is quite 
clear from the Figure 1 that the residues at the peak positions of the imaginary part are 
indeed positive. It is this coarse-grain picture that leads us to believe that all is well with the 
model and the predicted numerical values for the meson masses and decay constants are in 
agreement with the experimental values to within 10 — 20%. However, as is clear from Figure 
3, the imaginary part of the function (the spectral density) acquires negative values. Also, 
it is clear from the figure that we encounter negative values in a somewhat periodic manner 
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Figure 3: Im[f(x)] plotted in different intervals 



and that the magnitude of the negative values attained diminishes as we go to larger values 
of the argument. This feature is also not hard to expect and understand. The "somewhat 
periodic" appearance of the negative values is due to the oscillatory behaviour of Bessel 
functions (or their combination) while from the behaviour of the Bessel functions, it is very 
evident that the amplitude keeps on decreasing for larger and larger values of the argument. 

It is now straightforward to convince ourselves that the spectral density obtained within 
these models is not positive semi-definite. However, as we have discussed earlier, the spec- 
tral density is a positive quantity and has been very well measured in experiments, in full 
conformity with our theoretical expectations. This is the most important observation and 
result of this study. It may be worthwhile to briefly comment on the results/findings of the 
models [cOE] in the light of this observation. As has been mentioned above, if we just content 
ourselves with the positions of the resonances (meson masses) and evaluate the residues at 
these positions (decay constants), we'll obtain positive values for them. This is essentially 
what is done in [SJ E] • It is only when we take a closer look at the function under consid- 
eration that we are led to the observation that the imaginary part of the function, which is 
supposed to be positive always, acquires negative values. 

The above observation regarding the spectral density becoming negative is a serious issue of 
concern. As was mentioned, the spectral density is directly related to the hadronic r decay 
width by the optical theorem. A negative spectral density will simply mean that the decay 
width is becoming negative - something that is clearly unphysical and of course unobserved. 
Translated back, this observation has to say something very important about the model 
itself. In this form, the model does not lead to physical predictions and therefore can not 
be trusted. Similar conclusions can be reached at from the data on e + e~ annihilating into 
hadrons or the hadronic contributions to (g — 2) M . In each case the spectral density is related 
to a positive physical observable like cross-section or decay rate. 

We make a brief attempt to discuss the possible origin of such a problem in these models. 
Recall that to capture the essential features of chiral symmetry breaking, an ad-hoc and 
artificial infra-red boundary was introduced in the theory and certain specific boundary con- 
ditions specified on it. This approach is completely phenomenological and though, intuitively 
may seem well motivated and correct, by itself, does not guarantee that the results will be 
unitary and physical. As was pointed out initially, in a complete microscopic description, 
the geometry of the space-time should be appropriately modified and this should be consis- 
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tently done so that away from the AdS boundary (z — 0), the model incorporates the correct 
pattern for chiral dynamics. However, in the case at hand, this was avoided by invoking the 
artificial IR boundary. In our opinion, this itself is the root cause of the problem. The reason 
is as follows. Both the masses and the decay constants are obtained by essentially looking at 
the constant "b" that appears in the solution to the equation of motion for the gauge field. 
This constant is fixed in the present scenario by requiring the derivative of the solution to 
vanish at the artificial IR boundary, thus yielding a ratio of two Bessel functions. It is this 
combination of Bessel functions that leads to the trouble. We may be led to speculate that 
if, instead of the approximate form for the relative constant h in Eq([7jl. we had used the full 
expression, we would have bypassed the problem. It is again easy and straightforward to 
convince ourselves that this is not the case as the new form is also a combination of some 
other Bessel functions. If however, we can model the chiral symmetry breaking pattern by 
continuously changing the geometry in the bulk, this problem can possibly be avoided. We 
would like to point out that this is a common problem in all the models that invoke an 
IR boundary condition to model the chiral symmetry breaking. We restricted ourselves to 
the vector sector of the theory but the same arguments and conclusions apply to the axial 
vector sector as well. One can also check that the Weinberg sum rules are not satisfied in 
these models and there is no explanation of the (approximate) p meson dominance that is 
observed in nature. 

We conclude by summarizing our main observation and some of its consequences. We have 
investigated the profile of the spectral density in the recently proposed holographic models 
of QCD or chiral symmetry breaking in five dimensions [Sill]- The models seem to be quite 
robust and naively taken, seem to be predictive to within 10% accuracy. However, a closer 
look reveals that the spectral density, extracted from the two point correlation function, keeps 
acquiring negative values. This is in contrast to the positive behaviour of the spectral density 
expected from very general field theory arguments. Moreover, the vector and the axial- vector 
spectral densities are directly measured, for example, in the corresponding hadronic r decay 
modes. A negative spectral density implies a negative decay rate, in clear violation with 
unitarity and optical theorem and also with the observed data. This simply implies that the 
proposed models, though seem remarkably predictive, do not satisfy some of the basic field 
theoretic requirements and therefore can not be trusted. Also, in this form they can not be 
seriously taken to be models describing the physical hadron spectrum. Similar arguments 
will hold for any other model that violates the positivity condition for the spectral density 
and/or is in conflict with generic field theoretic expectations. As pointed out, the root cause 
in the present case is the way chiral symmetry breaking has been modeled by introducing 
an ad-hoc IR boundary. 2 If on the other hand, this is done in a more consistent manner by 
suitably modifying the bulk geometry in a continuous fashion, there is hope to get a dual 
model of hadrons which avoids the above mentioned problem. 
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need not apply to the works listed in J2J. 
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